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ABSTRACT 

We propose that the shape of the upper-end X-ray luminosity function (XLF) observed in elliptical 
galaxies for point sources carries valuable information about the black-hole (BH) mass spectrum among 
old X-ray transients formed in the galaxies. Here we present the line of arguments and analysis that 
support this connection and the methodology for deriving the BH mass spectrum slope from the 
observed XLF slope. We show that this underlying BH mass spectrum is modified by a weighting 
factor that is related to the transient duty cycle and it generally depends on the host-galaxy age, 
the BH mass and XRB donor type (main-sequence, red-giant, or white-dwarf donors). We find that 
the observed XLF is dominated by transient BH systems in outburst (a prediction possibly testable 
by future observations), but that the assumption of a constant duty cycle for all systems leads to 
results inconsistent with current observations. We also find that the derived BH mass slope depends 
on the strength of angular momentum loss due to magnetic braking for main-sequence donors. More 
specifically, we find that, for "standard" magnetic braking, BH XRBs with red-giant donors dominate 
the upper-end XLF; for weaker magnetic braking prescriptions main-sequence donors are found to be 
dominant. The methodology presented here can be used in the future as our understanding of the 
transient duty and its dependence on binary and mass-transfer properties improves. Under certain 
assumptions for this dependence we derive a differential BH mass spectrum slope of ~ 2.5; an upper 
BH mass cut-off at ~ 20 M Q is needed to understand the very brightest of the BH XRBs in elliptical 
galaxies. We also show that our quantitative results are robust against expected variations by factors 
of a few of the outburst peak X-ray luminosities. We expect that our analysis will eventually help to 
constrain binary population synthesis models and the adopted relations between black holes and the 
masses of their progenitors. 

Subject headings: galaxies: elliptical - methods: statistical - X-rays: binaries 



1. INTRODUCTION 

Chandra has revolutionized the study of point X- 
ray sources in the nearby Universe. The majority of 
these are interpreted to be X-ray binaries (XRBs; for a 
general review on Chandra performance see Weisskopf 
et al. 2003, for extragalactic X-ray binaries see, e.g., 
iKim fc Fabbianol 12004 Uordan et all 120041) . Elliptical 
galaxies out to the Virgo cluster have now been stud- 
ied and have surprised us with the large number (typi- 
cally ~ 100 per galaxy) of detectable point X-ray sources 
down to X-ray luminosities of about 10 37 ergs _1 . Two 
main population characteristics have attracted consider- 
able attention so far: (i) the X-ray luminosity function 
(hereafter XLF) that may or may not exhibit a break 
at about 4 — 5 x 10 38 ergs~ 1 (for a recent update see 
IKim fc Fabbianol |2004|) ; (ii) the high fraction of sources 
coincident with identified globular clusters (GCs) in el- 
lipticals. 

The shape of the XLF has been debated since the first 
observ ations of ellipticals were reported. iSarazin et alJ 
(2000) identified a shape that required two power laws 
with a "break" or a "k nee" at ~ 3.2 x 10 38 ergss _1 . 
IKim fc Fabbianol l)2003j) argued that the break may re- 
sult from biases affecting the detection threshold of the 
data. In the following few years longer exposures be- 
came possible and more and more ellipticals were added 
in the observed sample with low enough sensitivity (see, 
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e.g. iGilfanovl 120041 ) The current situation is probably 
best summarized in IKim fc Fabbianol l|2004fl . They an- 
alyzed a large sample of elliptical galaxies with varying 
sizes of point source samples, and they concluded that: 
although XLFs of individual galaxies do not require a 
broken-power-law fit, the combined sample of sources 
from all the galaxies considered shows a statistically sig- 
nificant requirement for two power laws and a break at 
5±1.6x 10 38 ergs s _1 . They found the best-fit slope of the 
lower end of the differential XLF to be ad = 1.8±0.2 and 
the best-fit slope of the upper end to be ad = 2.8±0.6. It 
is important to note that the break location is consistent 
with the Eddington luminosity for a 1.9±0.6M Q neutron 
star (NS) accreting helium-rich material (for hydrogen 
rich donor this value is as large as 3.2 ± 1M W ). In what 
follow s we consider the results of the IKim fc Fabbianol 
(2004) study as representing our current observational 
understanding of the XLF in ellipticals. We address the 
question of the interpretation of this understanding and 
what it implies about the properties of the sources con- 
tributing to the observed XLFs. 

Large fractions (20% - 70%) of the point sources in 
ellipticals have been r eported to be associ ated with glob- 
ular clusters (see, e.g. . ISarazin et aH2003l and references 
therein). These high fractions have led to the suggestion 
that all point X-ray sources seen currently in ellipticals 
have been formed through stellar interactions and that 
sources that are not associated with GCs have originated 
in GCs and have either been ejected or the parent GCs 
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have been destroyed by the galaxian tidal forces. As 
much interesting as this suggestion is, it raises the ques- 
tion: why would the field stellar population of ellipticals 
not lead to XRB formation as it has occurred in the 
Milky Way, for example? One could speculate that the 
field population is just too old and the XRBs that were 
formed at some point have completed their X-ray emit- 
ting life. However, such a speculation is inconsistent with 
the expectation that XRBs with low mass donors can live 
for several Gyrs as the donors lose mass and the binaries 
enter a transient phase. Such systems would become de- 
tectable as bright X-ray sources during the disk outbursts 
ijPiro fc Bildstenl 12002ft . Moreover, it has recently been 
pointed out that the high rate of source coincidence with 
GCs actually appears consistent with some of the sources 
having been formed in the field ijJuettl 12005}) . although 
the result may sensitively depend on the definition of 
GC concentration in ellipticals. Most importantly for 
the present study, of the bright sources above the XLF 
break only a very small fraction is associated with GCs 
(e.g., only one source in the Virgo cluster sample as re- 
ported bv I Jordan et al.ll2004D. 

In a recent Letter iBildsten & Delovel l)2004 hereafter 
DB) have suggested an explanation that couples the 
two population characteristics: the XLF shape and the 
source coincidence with GCs. They identify the point X- 
ray sources as ultra-compact binaries (UCBs) that form 
predominantly in GCs. They are neutron stars accret- 
ing from low-mass He or C/O white dwarfs and they 
contribute to the ellipticals XLF early in their lifetime 
when they are still bright. Their association with GCs is 
important in replenishing the population through tidal 
interactions and allowing a significant number of sources 
in this bright phase, even though there is no ongoing 
star formation in ellipticals nor in GCs. DB estimate the 
rate of ultracompact binary formation to be consistent 
with the number sources observed and conclude that the 
model XLF slope is i n agreement with the sl ope below 
the break as derived bv lKim fc Fabbianol (2004). The as- 
sociation of the XLF break with the NS Eddington limit 
for He-rich accretion is also consistent with this inter- 
pretation. However, the upper end of the XLF (at lumi- 
nosities in excess of the break location) are not easy to 
interpret. Deloye & Bildsten suggest that some of the NS 
ultra-compact sources can reach super-Eddington lumi- 
nosities. However luminosities in excess of 10 39 ergs _1 
are very difficult to explain with NS accretors. There- 
fore the origin of the upper-end slope is not naturally 
connected to NS UCBs formed in GCs. 

In this paper we address the question of the upper- 
end XLF slope and its origin. We consider the pre- 
viously made suggestion ijSarazin et alJ 12000(1 that the 
XLF above the break at 5 x 10 38 erg s" 1 is populated by 
XRBs with black hole (BH) accretors. Given that GCs 
are not ex pected to harbor a si gnificant number of BH- 
XRBs fsee lKaloeera et al"ll2004 and references therein), 
we suggest that the vast majority of these BH-XRBs are 
part of the galactic-field stellar population in ellipticals. 
As we will show most of donors in these binaries are of 
low-enough mass that the XRBs are expected to be tran- 
sient and therefore they populate the XLF only during 
disk outbursts when they typically emit at the Edding- 
ton luminosity for their BH mass. We further suggest 
that the slope of the upper XLF is a footprint of the BH 



mass spectrum in the BH XRBs under consideration. We 
present analytical derivations that demonstrate this link 
and we develop a method that allows us to infer the un- 
derlying BH mass spe ctrum consistent with t he current 
upper-end XLF slope l)Kim fc Fabbiano| |2004L We also 
show that given the current observations it is possible 
to constrain the strength of magnetic braking acting in 
these XRBs, the type of BH donors, as well as the tran- 
sient duty cycle to some extent. We also examine the 
quantitative robustness of our results against variations 
of some basic assumptions. This analysis is presented in 
§ 2 and 3. We conclude with a discussion of our results 
and possible connections to population synthesis calcu- 
lations (§4). 

2. BLACK HOLE X-RAY BINARIES IN ELLIPTICALS 

We consider XRBs that could possibly populate the 
part of the observed XLF above the reported break at 
4 — 6 x 10 38 erg s _1 , and therefore we focus on BH ac- 
cretors (masses in excess of 2 — 3M©). Given the current 
estimates for the ages of stellar populations in ellipti- 
cals (in their majority 8 to 12 Gyr, althoug h some es- 
timates are slightly sh orter than 5 Gyr; see lRvden et alJ 
l2001tlTemi et al.ll2005fl . we expect that donor masses are 
lower than ~ 1 — 1.5 M©. Given these mass ratios and 
the properties of similar observed systems in the Milky 
Way (i.e., soft X-ray transients), these BH XRBs are ex- 
pected to be transient X-ray sources (see McClintock & 
Remillard 2005 for a review of BH X-ray binaries in the 
Milky Way) , where mass transfer is driven by the Roche- 
lobe filling donor. Given the above upper limit on the 
donor mass for ellipticals, we expect that there will be 
three different types of low-mass donors: (i) Main Se- 
quence (MS) stars, (ii) Evolved or Red Giant Branch 
(RG) stars, and (iii) White Dwarf (WD) donors. Each 
of these sub-populations of BH XRBs will have different 
typical mass-transfer rates and binary property distribu- 
tions, and therefore we examine them separately in our 
analysis that follows. 

2.1. Transient X-Ray Sources 

We adopt the current understanding for the origin of 
transient b ehavior in XRBs (for a recent review, see 
King 2005). To identify transient systems in our mod- 
eling we consider the typical mass-transfer (MT) rate 
associated with each type of XRB donor (Mi) and com- 
pare it to the critical MT rates for transient behavior 
(-^fcrit): if Mi < M cr it, then the accretion disk is ex- 
pected to be thermally unstable and the binary system 
is assumed to be a transient X-ray source. The value 
of this critical MT rate for the disk instability is not 
precisely known and its functional dependence on disk 
and binary properties are subject to uncertainties as- 
sociated with our current theoretical understanding of 
the disk instability. However, both the qualitative con- 
cept of the existence of a critical rate for the instability 
to set in and its quantitative estimates by recent stud- 
ies appear to be in good agreement with the behavior 
of Galactic X-ray transients. Therefore, we adopt the 
current understanding and quantitative estimates. More 
specifically, for h ydrogen-rich donors , we adopt the M cr it 
value derived bv lDubus et al.l (^999), and for helium or 
carbon-oxygen donors, we adopt the value derived by 
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M enou et aT] (J2002). The effects of quantitative devia- 
tions from the adopted expressions are discussed in what 
follows. 

To account for the contribution of transient sources 
in the XLF among any persistent sources, assumptions 
about the XRB luminosity at outburst and the transient 
duty cycle need to be made. 

When a XRB is identified as transient, we assume that 
during the disk outburst the X-ray luminosity Lx is equal 
to the Eddington luminosity L B dd associated to the BH 
accretor: 

L E dd = =5xl(T' ergss L , (1) 

K K 

where mbh is the accretor mass in M@ and k is the 
opacity of the accreting material in cm 2 g _1 . We adopt 
electron scattering opacities equal to 0.32 and 0.19 for 
hydrogen and helium or carbon-oxygen rich material, re- 
spectively. 

We note that of the 15 confirmed transient BH XRBs in 
our Galaxy, 3 appear to reach possibly supe r-Eddington 
luminosities llMcGlintock fc Remillardll2005D at outburst 
(although distance estimate uncertainties cannot be ig- 
nored). Two of them, V4641 Sgr and 4U 1543-47, have 
early-type donors. Such donors are not present in ellip- 
tical galaxies with population ages of ~ 5 — 10 Gyr. The 
third one, GRS 1915+105, has a low-mass giant donor 
and an orbital period of 33 days. However its X-ray lumi- 
nosity at outburst just barely exceeds its LEdd, by 40% 
only. Given distance uncertainties associated with such 
an estimate, we conclude that we can neglect the pos- 
sibility of super-Eddington luminosities during outburst 
in our XLF modeling. On the other hand outburst peak 
luminosities cover a significant range at sub-Eddington 
values. During primary 1 outbursts peak Lx values can 
be lower than L B dd by factors of a few (McClintock & 
Remillard 2005, private communication). As part of our 
analysis we examine the effect of such variations on the 
methodology and conclusions presented here (see § 3.1.1). 
From an observationally point of view it has been shown 
(Zezas et al. 2004 and Zezas 2005, private communi- 
cation) that variability in X-ray fluxes (and hence lu- 
minosities) by factors of a few (typical among accreting 
sources and detected with Chandra observations at dif- 
ferent epochs) do not alter the XLF slopes as measured 
for nearby galaxies within the current errors. Therefore 
observationally the reported XLF slopes appear to be 
robust. Consequently we can use them to learn about 
the underlying XRB population with considerable confi- 
dence. 

At present there are no strong constraints on the duty 
cycles either from observations or from theoretical con- 
siderations. Among known Galactic X-ray transients, 
typical duty cycles of a few % is favored for hydrogen 
donors l|Tanaka fc Shibazak3 ll99^1 . To our knowledge, 
there are no data on duty cycles for transients with a WD 
companion. In what follows we investigate how plausible 
duty cycle assumptions affect the upper-end XLF shape. 

We use the term "primary" to distinguish from "follow-up" 
outbursts that are occasionally observed very soon after primary 
ones with peak luminosities orders of magnitude below the Edding- 
ton limit (Remillard & McClintock 2005, private communication). 
Such small outbursts do not reflect an extremely short duty cycle 
and do not contribute to the high-end XLF of interest here. 



In particular, we consider two specific cases: one of con- 
stant duty cycle equal to n = 0.01; another of a variable 
(dependent on MT rates) duty cycle equal to 

l-M f^-V (2, 



M, 



crit 



where 5 = 1 is assumed. The first of these two cases 
corresponds to the standard assumption of a constant 
duty cycle often made in the literature. The second case 
is motivated primarily by our plan to examine how one 
example form of a MT-dependent duty cycle affects our 
analysis and results. Admittedly the specific choice of 
the dependence on M CT it shown above is not solidly mo- 
tivated, given all the uncertainties of the outburst mech- 
anism. However it implies a correlation of the duty cycle 
with how strong a transient the system is: the further 
away from the critical MT rate, the smaller the duty cy- 
cle. We stress that in most of our analysis we adopt this 
form with 5 = 1 as a plausible example and throughout 
the paper we contrast the results to those obtained with 
a constant duty cycle. Furthermore in § 13.1.11 we exam- 
ine the sensitivity of our results for main-sequence donors 
on the choice of the duty-cycle dependence on MT ex- 
tensively: we adopt 5 > 1 in eq.© and we also introduce 
yet one other example of a MT-dependent duty cycle: 

.-»(&)' 

for which we examine various values of S. Once again 
there is no solid theoretical motivation for this latter 
functional choice. However, it provides us with a better 
understanding of how sensitive our results are to the de- 
tails of the possible duty-cycle dependence on MT prop- 
erties. 

2.2. Main Sequence Donors 

Mass transfer in BH XRBs with hydrogen-rich, low- 
mass MS donors is expected to be driven by angular mo- 
mentum losses due to magnetic braking (MB) and gravi- 
tational radiation (GR) . In the case of conservative mass 
transfer l|Verbuntll 19931) the angular momentum loss rate 
are connected to the MT rate as follows: 



Jorb 



M d V 6 



n 
2 



Al 



u / - (4) 

where n is the radius-mass exponent for the donor. For 
a low-mass MS star: 

ra — md, (5) 
where r d = R^/Rq and rn d = M^/Mq are the donor 
stellar radius and mass in solar units. Therefore n = 
d InRd/d InMd is equal to ~ 1 for MS donors. 

According to general relativity the rate of angular mo- 
mentum loss due to GR is given by: 



Jorb 



32G 3 M BH Md(M B H + M d ) 



5c 5 



,4 4 



= -2.6 x 10- 17 mBHmd(w 4 BH+Wd) s-\ (6) 

a 4 

where a is the orbital semi-major axis in un its of solar 
radiu s. For a mass ratio q = M^/Mbh < 0.8 ( Pa czvnskl 
GUZI and using the mass-radius relation eq. (JSJ): 



1 2 /3/ , U/3 

— ml (m BH + m d ) 1 . 



(7) 
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We consider two derivations of the angular momentum 
loss rate due to magnetic braking: (i) the Skumanich law 
based on the empirical relation for sl owly rotating stars 
adopted from (fRappaport et alTl 983) (RVJ) , and (ii) the 
revised law based on X-ray observations of fas ter rotating 
dwarfs adopted from l)Ivanova fc Ta am 2003) (IT): 

iT= -3-8 x W- w M d R^(R d /R G ) 2 n 3 dyncm (8) 
JZ = -6 x 10 30 (i? d /i? Q ) 4 (^p) dyncm, (9) 

where [s _1 ] is the stellar angular velocity which is equal 
to the binary orbital velocity assuming the star is in full 
synchronization with the binary orbit, fi g = 5 x 10~ 6 
s _1 is the Sun's angular velocity, and £! x = 10il Q . Using 
Kepler's law the above are re-written as: 

SL = - 7 .2 x 10- 15 m * (mBH + md)2 s- 1 , (10) 
J orb m B Hfl 5 

ffi = -2.7 X 10- 17 W ' (mBH + ™f s- 1 . (11) 

In XRBs with BH masses > 3M© and MS donor masses 
< 1.0M Q , it is J gr < J^ J and J gr > J™ h . The lifetime 
of the MS-BH XRBs is much longer in the latter case. 

As mentioned earlier in this study we adopt the deriva- 
tion of the critical MT rate below which the accretion 
disk becomes unstable for irradiated disks presented by 
Dubus et al. 1999: 

M crlt = -1.5 x lO 15 ™^ 4 ( I ^) 21 gs- 1 . (12) 

Here -Rdisk is the radius of the accretion disk. We note 
that the exact value of this critical rate is subject to 
uncertainties associated with our limited understanding 
of the disk instability, but we adopt the above expression 
as indicative of the process and we continue with our 
analysis. 

From eq.|g|), Un|) and O it can be shown numerically, 
that for the RVJ MB law and for low-mass donors, there 
is a BH mass Mpt of ~ 5 M Q that separates BH-MS 
systems into persistent (Mbh < Mpx), and transient 
(Mbh > Mpt)- From more detailed binary evolution- 
ary calculatio ns using the stellar evol ution and MT code 
described in ijlvanova fc Taarr] 120041) . we find that this 
boundary is about 10 M Q (see Fig.^for details). For BHs 
less massive than this critical mass, the XRBs are persis- 
tent as long as the donor masses are higher than about 
0.3 Mq. In these persistent sources the MT rates driven 
by the RVJ type of MB turn out to be 0.01 - 0.25 of the 
black holes's Eddington rate. As a result, the persistent 
X-ray luminosity for these systems is < 10 38 ergs _1 , i.e., 
below the bright Lx range we consider here. Therefore 
we conclude that the persistent BH-MS binaries driven 
by the RVJ type of MB cannot contribute significantly 
to the upper-end XLFs of ellipticals. 

Next we examine whether the transient phases as- 
sociated with BH-MS binaries and the RVJ MB law 
are important when they reach X-ray luminosities com- 
parable to the Eddington limit. For Mbh > 10 M , 
the outburst luminosity is expected to be in excess of 
~ 1.5 x 10 39 ergs _1 . However, this lower limit is compa- 
rable to the highest luminosity seen currently in XLFs of 
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Fig. 1. — Evolution of BH-MS binaries with the RVJ magnetic 
braking prescription, for BHs of 3 and 10 Mq; the initial MS com- 
panion mass is 1 Mq. Shown are the MT rate M (solid line), the 
critical MT rate Af cr ; t (dotted line) and the Eddington MT rate 
Mgdd (dashed line), all rates are in Mq per yr. The dash-dotted 
line shows the donor mass evolution. 



ellipticals (Ki m &: Fabbianol l2004|) . and therefore these 
systems cannot contribute significantly to the observed 
XLFs. The last possibility is outbursts from transient 
BH-MS with Mbh < 10 Mq and donors less massive 
than ~ 0.3 M . Such low mass donors are out of ther- 
mal equilibrium and significantly expanded (~ 3x) com- 
pared to an undisturbed MS star of the same mass. In 
the case of the MT dependent duty-cycle ij is about a 
few %. We note, however, that applicability of MB for 
these stars is very questionable, as it is generally accepted 
that MB does not operate in fully convective stars, which 
are found to be less massive than 0.35 M Q for undis- 
turbed stars. In principle, however, "eroded", out-of- 
thcrmal-cquilibrium low-mass MS donors like the ones in 
XRBs do not necessarily become fully convective at the 
same critical mass as stars with no prior MT evolution. 
For this reason we have used detailed MT calculations 
(Fig.^l with a stellar-evolution code to examine this be- 
havior further. We find that in the donor-mass range 
0.15 — O.3M the radiative core is extremely small, even 
for these "eroded" stars, and therefore applying angular 
momentum loss due to MB is not reasonable. Instead 
MB activity and hence mass transfer is expected to be 
interrupted until eventually GR drives Roche-lobe over- 
flow much later. 

Based on the above line of arguments we conclude that 
BH-MS binaries evolving according to the RVJ MB law 
are not expected to contribute significantly to the high- 
end XLFs of ellipticals. 

In the case of the IT MB prescription, BH-MS systems 
are transient for all BHs masses A/bh > 3M and for all 
low-mass MS donors. The reason is that the IT MB is 
weaker; consequently the donors are mildly out of ther- 
mal equilibrium and the mass transfer rates are lower 
compared to the RVJ MB case. Using again detailed 
MT evolutionary simulations we find that M/M cr ;t ~ 
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Fig. 2. — Evolution of BH-MS binaries with the IT magnetic 
braking prescription, for BHs of 3 and 10 Mq; the initial MS com- 
panion mass is 1 Mq. Shown are the MT rate M (solid line), the 
critical MT rate M cr ; t (dotted line) and the Eddington MT rate 
-Mgdd (dashed line), all rates are in Mq per yr. The dash-dotted 
line shows the donor mass evolution. 



0.25 ± 0.15 (see Fig. Therefore, in the case of the 
MT dependent duty-cycle r\ is again about a few %. In 
both cases (RVJ MB and IT MB), the value of the duty 
cycle is consistent with observations for BHs of different 
masses, though the transiency occurs at very different 
donor masses. 

We conclude that, regardless of the specific MB law, 
it is rather unlikely that persistent sources with a BH 
accretor and a hydrogen-rich, low-mass MS donor pop- 
ulate at any significant fraction the upper-end XLF of 
ellipticals; only transient BH-MS sources driven by the 
IT MB law can populate this X-ray luminosity range. 

2.3. Red Giant Donors 

For orbital periods more than about a day, MT occurs 
when the low-mass donor is a subgiant or a giant. The 
driving force is the nuclear expansion of the donor, and 
a simple analytic prescription for the MT is (Webbink, 
Rappaport, Savonije 1983; Ritter 1999; see also King 
2005): 

m 1.47 



Mrg = -3.4 x 10 



15 a 1A 



gs 



(13) 



(TOBH + TO d ) - 465 

It has been shown Cbv lKing et aHll997llKinell200ffl that 
such wider XRBs are transient, regardless of the BH mass 
(original derivations were based on a somewhat different 
expression for the critical MT rate for transient behavior, 
but still quite similar to eq. |12j . 

2.4. White Dwarf Donors 

A typical WD ma ss-radius relation is (see, e.g., 
iRapnanort et al]ll987|) : 



r d = 0.0128m 



-1/3 



(14) 



Using an appro ximation for the Roche Lobe radius (from 
IP acz vhskl 1 1 9 7 1[) and assuming that the mass of WD is 



much smaller than a BH mass, we can show that 

-1/3 



0.0278 m; 



1/3 / ra d 



(15) 



We consider again conservative mass transfer J3J but 
without any MB losses, adopting n = —1/3 and assum- 
ing that m-d <C m,BH- Then 

32G 3 M 2 „M 2 



BHi " d = -7 9 X lQl6'"BH'"d - 

3.3c 5 A 4 r 4 
We substitute here eq. \ib\ and then have 

'bh M ©y r 



-3 4" 

md = — 2 x 10 m^mi 



(16) 



(17) 



In what follow s we adopt the critic al MT rate for He- 
rich donors from lMenou et aD lj2002fl . but we note that 
the expression is subject to quantitative uncertainties as- 
sociated with the current understanding of the disk in- 
stability: 



Mr. 



-5.9 x lO 1 ^ 87 (^f 



gs 



(18) 



For a large mass ratio (/bh = WBH/wd the Roche lobe of 
the accretor is trx — 0.7a and 



fdisk ^ 2/3r RL 



0.013m d 2 ^rn^^ 



m crit = -1.7 x 10- 12 mJ 174 M^r" 1 
BH-WD binaries will be transient if 
m d 



1.2 x 10 9 m^- 4 mj| H < 1 



(19) 
(20) 

(21) 



"Icrit 

Therefore the maximum donor mass that leads to tran- 
sient behavior in BH-WD binaries is: 

-o.i 



mtr = 0.038 m 



bh 



(22) 



The time interval in Gyr needed for the WD donor 

i_ ]_/ 3 11/3 

mass to evolve from m d (Ti) to m d (T2) is (using 
eq.EJl: 



T 2 - T x : 



3 
22 



x 10" b m 



-2/3 

BH 



md(T 2 ) 



-11/3 _ 



md(Ti) 



-11/3 



(23) 



Here we assume that the mass of the BH is constant, 
since wbh 3> mj. Consequently and using eq. (|22|l we 
can find that the time a BH-WD system spends in the 



persistent state is t pe 



20 x 10 b m 



-0.3 

BH 



yr, i.e., it very 



weakly depends on the accretor mass (the dependence on 
the initial donor mass is negligible, below that 1%). We 
note that through this persistent phase the MT rate will 
be comparable or higher to the Eddington limit only for 
a very short time, ^Edd — 2 x 10 6 m.g H ' 9 yr 2 . 

The evolution of BH-WD systems with C/O WD com- 
panions is rather similar. The critical MT rate (using 
Menou et al. 2002) is 



m crit = -9.4 x 10 



-13, 



- 1 - 47 Moyr" 1 



(24) 



2 Although MT is non-conservative during the super-Eddington 
accretion, and eq. J1J formally should not be applied, this result is 
well consistent with the detailed calculations that take into account 
non-conservative MT 
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and 

m tr = 0.03 mgH 1 • (25) 

The time that a BH-WD system with a C/O rich donor 
spends in the persistent state is also not very long, t pcrs ~ 
50 x 10 6 m-° 3 yr. 

We conclude that BH-WD binaries that contribute to 
the current upper-end XLFs of ellipticals are expected to 
be transient sources 3 

3. MASS SPECTRUM WEIGHTING FACTOR AND 
TRANSIENT DUTY CYCLE 

In the previous section we have shown that the upper- 
end XLF of ellipticals is dominated by transient BH 
XRBs possibly with a variety of donors: MS (for the 
case of the IT MB prescription) and RG stars, and WD 
donors with masses lower than ~ 0.035 Mq. All these 
systems contribute to the XLF only when in outburst, 
when their Lx — ^Edd °c Mbh- Consequently the slope 
of the upper-end XLF can serve as a footprint of the 
BH mass distribution of accretors in the contributing BH 
XRBs. These contributing BH XRBs are just a sub-set 
(those in outburst) of the true population of BH XRBs in 
ellipticals determined by the duty cycle of BH transients 
binaries. For the general case of a transient duty cycle 
that is dependent on the BH accretor mass (and possi- 
bly other quantities), the differential XLF n(L) ^ s and 
the underlying BH mass distribution in XRBs n(m)BH 
are connected by: 



n(L x )obs = n(m BH ) x W(m B n), 



(26) 



where VF(m BB ) is a weighting factor related to the de- 
pendence of the transient duty cycle on m BB - 

The observed slope of the differential upper-end XLF 
is ad = 2.8 ± 0.6: n(Lx)obs oc L x ad (the slope of the 
cumulative upper-end XLF reported by Kim & Fabbiano 
2004 is a c = 1.8 ± 0.6). Assuming that n(m B H) oc m B {| 
and W(mBB\) oc m^, the slope characterizing the un- 
derlying BH mass distribution in XRBs is: 



P = a d - 7- 



(27) 



For the standard assumption of a constant duty cycle, 
/3 = ay = 2.8±0.6. In the following subsections we 
derive W(mBn) and 7 for all three types of donors in 
one example case of a duty cycle dependent on the binary 
and MT properties (see, e.g., eq0). 

3.1. Main Sequence Donors 

In what follows we estimate the typical duty cycle for 
BH-MS transients averaged over the possible distribution 
of donor masses. In the case of the IT MB perscription, 
the angular momentum loss rate due to GR is compara- 
ble or even more important than MB for all BH masses 
above 3 M Q and donor masses < 1.0 - 1.2M©. The MT 
timescale during the transient phase is longer than the 
donor's thermal timescale when the donor is > 0.25 M Q . 
For this range the donor is in thermal equilibrium and the 

3 This would not be true if BH-WD binaries continuously formed, 
but this is not possible in the galactic field of ellipticals and is 
not even expected in globular clusters, since BHs tend to dynam- 
ically separate from the rest of the cluster and eject one a nother 
(Kuikarni ct al. 1993; Sijjurdsson & Hcrnauist 1993; Wattcrs ct al. 
120001) . 



approximation for the mass-radius dependence eq. JSJ 
can be used. Using eq. JIJl, @) and ijT21) , and the fittin g 
formula for the Roche lobe radius from lEggletori ((19831) . 
we find: 

™ 05 i w BH(gBH /3 Ml + £bh) + 0-6) 21 , , 
m crit ' TO2(l + gB H)(4/3-l/ gB H) ' 1 ' 

where qbh = ^BYi/md is the mass ratio. For g B H ^ 1 
we have 



0.0Um B j l m d 



(29) 



At donor masses smaller than 0.25 M Q , the donor is 
out of the thermal equlibirum and its radius is about 
twice bigger than predicted by eq. We then find: 



m cri t 



0.0004m B -£m d 



(30) 



We note that the split into the two expressions above is 
a rough, but useful approximation. 

As discussed in §2.2, for IT MB, a BH-MS system 
is transient throughout the MT phase. The MT rates 
are well below the Eddington limit for the BH mass and 
therefore we assumed that MT is fully conservative; i.e., 
Md + A/eH = Mtat is constant with time. In what follows 
we assume a flat current mass distribution for donors 
(dN/dnid — const). We are guided in this choice by re- 
sults from binary population synthesis calculations (with 
the Star Track code; Belczynski et al. 2002 and 2005; 
Belczynski 2005, private communication). We integrate 
eq. [25]for from ~ 0.25tolM Q and using eq. we 
find that at present the probability that a system with 
a BH accretor of wbh is in outburst and therefore con- 
tributes to the upper-end XLF is: 



W(m B n) = 



It 



™to ON 
0.25 "dm d 



f "iTO ON 
Jo.25 a^J dm d 



0.05m B H, 



(31) 



where mxo is the turn-off MS mass for the elliptical 
galaxy in solar units. We note that this result is valid 
only for large mass ratios g BB >> 1. The contribution 
of BH-MS system when donors have masses < 0.25A/ Q 
(systems where the donor is out of the thermal equilib- 
rium) is less significant. 

It is also important to note here that for MS donors the 
factor W does not appear to depend on time (i.e., the age 
of the elliptical galaxy). Such a time dependence would 
enter in relation to the value of the maximum donor mass 
(turn-off mass for the host galaxy). However we find that 
W(mBn) is a very weak function of mxOi an d therefore 
it is not sensitive to the elliptical age. 

3.1.1. Monte Carlo Simulations 

In principle, prolonged mass accretion onto the BHs 
can affect their mass spectrum. Since this effect can- 
not be included analytically, we have examined it quan- 
titatively using simple Monte Carlo simulations. We set 
up the simulations assuming a flat BH-MS birth (MT 
onset) rate and a flat mass distribution for donors at 
the onset of the MT phase, without any restrictions on 
the mass ratio g B H- Donor masses at birth were var- 
ied in from 0.1 Mq to the current mxo a t an ellipti- 
cal age of 10 Gyr, assumed as a standard value. Each 
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Fig. 3. — The upper panel shows model XLFs for BH-MS 
binaries. Lines shows results for initial /3o = 1.7,2.3,2.9 at an 
age of 10 Gyr (dotted, solid and dashed lines; respectively). The 
dash-dotted line is for a 15 Mq BH mass cut-off with /3o = 2.3. 
Thick solid lines corresponds to slopes a c = 0.8 and ct c = 2.3. The 
lower panel shows XLFs in observed early type galaxies, where 
it is shown that the observed range of slopes is also within the 
range cr c = 0.8 — 2.3. Data are taken from Kim & Fabbiana |l2004l) 
and the ages of t he ellipticals are from iRvden et alj I2001B and 
ITemi et all <2005l) . 




Fig. 4.— Model XLFs for BH-MS binaries for initial f3 = 2.3 at 
an age of 10 Gyr. The upper panel includes results for rj given by 
eq. JzJ for <5 = 1 and outburst peak luminosity deviating by a factor 
of 2 from the Eddington luminosities (the dash-dotted line), and for 
5 = 0.5,1,2 without any luminosity variation (dotted, solid, and 
dashed lines, respectively) . The bottom panel includes results for 
r\ given by eq. f° r S = 0.5, 1, 2 without any luminosity variation 
(dotted, solid, and dashed lines, respectively). Straight lines with 
various a c slopes are also shown. 



BH-MS binary was evolved to the current elliptical age 
using equations shown in the § 2.2. For the MT evo- 
lution we took into account both IT MB and GR, and 
for the Roche lobe radius we adopted the approximation 
by Eggleton (1983). A binary is removed from the MT 
population if the donor mass falls below 0.05 M Q . If the 
MT timescale is longer than the thermal timescale of the 
donor, the donor radius evolution is simply proportional 
to the mass lost due to MT. On the other hand, if the MT 
timescale is shorter than the donor's thermal timescale, 
the donor is out of the thermal equilibrium. In this case 
we modify the evolution of the donor radius using a pre- 
scription that is in acceptable agreement with our results 
from detail ed MT ca lculations with the stellar evolution: 
6r ~ Sm^J toth/™, where ttith = ^d/^TH is the MT 
rate driven on the donor's thermal timescale £th- Evo- 
lution of the transient systems follows the adopted duty 
cycle (eq.Jlj with 6=1). 

Based on the results of our Monte Carlo simulations 
we find that: (i) due to accretion the BH mass spectrum 
slope increases by about 0.2, i.e., (3 — (3q + 0.2, where 
0o is the BH mass slope at MT onset; (ii) the slope of 
the BH mass spectrum at the beginning of mass transfer 
best reproduces the observations with $ = 2.3 ±0.6 (see 
Fig. EI) . We also find that the relation between f3 and 0o 
is not sensitive to the age of the elliptical (as long as it 
is in the range of a few to several Gyr). 

Next we examine how the results in this section are 
affected by plausible variations of a number of possibly 
oversimplifying assumptions made so far. Although it is 
possible to re-derive the analytical expressions with dif- 
ferent assumptions, it requires repeating essentially the 
same analysis multiple times, something inappropriate 
for presentation here. Instead it is more efficient to ex- 



amine these effects using the Monte Carlo simulations for 
XLF slopes. 

For the tests that follow we adopt (3 = 2.3 and examine 
how the cumulative XLF slopes a c are affected. 

We first examine the effect of random variations by a 
factor of 2 of outburst peak X-ray luminosities of indi- 
vidual transient systems. The results are shown in Fig. 4 
(top panel): the solid line is our standard model where 
peak Lx are set equal to ^Edd (same as the solid curve 
in Fig.[j|]top panel) and the dash-dotted line is the result 
we obtain when random Lx variations are introduced. 
They are essentially indistinguishable, and therefore the 
derivation of a BH mass spectrum slope from the ob- 
served XLF slope appears highly robust. This is con- 
sistent with the findings of observational studies of such 
variations (Zezas et al. 2004). 

Next we examine the effect of 5 values in the case of 
the MT-dependent duty cycle in eq.El other than or 1. 
A set of curves for 8 =0.5, 1, and 2 are shown in Fig. 4 
(top panel) . Once again we conclude that the variations 
are essentially negligible to such quantitative changes. 
Instead the slope behavior seems to be dominated by 
the qualitative character of rj in eq.|21 the stronger the 
transient character, the smaller the duty cycle. 

Last we examine the effect of changing the functional 
dependence of the duty cycle normalizing it to the Ed- 
dington MT rate instead of the critical rate for transient 
behavior (eq. [2j). The results are shown in Fig. 4 (bot- 
tom panel) for 8 =0.5, 1, and 2 and they are compared 
to the our standard case of rj (from eq. [2] with 6=1). 
This is the only case where significant variations are ev- 
ident. More specifically, the resultant XLFs are steeper 
(absolute a c values are higher) and they are more depen- 
dent on the choice of 8. Given that the results are shown 
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for a fixed (3 value, it means that, for a specific observed 
XLF slope, this different form of the MT-dependent duty 
cycle would lead to the derivation of a flatter BH mass 
spectrum compared to that derived using eq. Q . We 
conclude that better understanding of the transient duty 
cycle and its dependence on binary or MT properties is 
important for obtaining reliable quantitative results in 
the future. 

3.2. Red Giant Donors 

In principle we should repeat the above estimate for 
BH-RG transients. However, from eq. 10, (|13|) and 
(|19[) we obtain: 



given BH mass: 



m C rit 



0.2a-°- 7 m}, 5 , 



(32) 



and it is evident that, in this case, the MT-dependent 
duty cycle 77 (see eq. 01 is independent of the BH mass. 
Therefore, for RG donors, 7 = 0. Instead it is signifi- 
cantly dependent on the RG donor mass. Since in ellip- 
ticals the typical mass for RG donors is about the same 
as the mass of the turn-off of MS stars, the duty cycle 
for RG donors depends on the turn-off mass, and hence 
on the age T [Gyr] of the elliptical. For a solar metal- 
licity and stars of 1M© < M < 1.5 M©, the approx- 
imate dependence of the turn-off mass with the age is 
m d ~ mxo ~ 2T -1 / 3 (here we used th e simplified evo- 
lutionary code from iHur lev et al1l2000|) . Using eq. (0), 
and integrating eq. (|32|) over binary separations (similar 
to our integrations for MS donors; see en. IHljl. we find: 



W(T) ~ 0.03 T 



-0.5 



(33) 



Here we assume a distribution of orbital separations for 
BH-RG binaries before MT starts that is flat in the log- 
arithm. The reason for this choice is that this is ap- 
propriate for the distribution of zero-age binaries and 
the shape is actually preserved through wind mass loss, 
common-envelope evolution, and asymmetric explosions 
with small kicks appropriate for black holes (Kalogera & 
Webbink 1998). 

It is interesting to note that for RG donors W is de- 
pendent on the galaxy age whereas for MS donors the 
dependence on the BH mass dominates. 

3.3. White Dwarf Donors 

During the transient stage, the WD mass can be writ- 
ten as a function of time T in Gyr (using eq. I23fl : 



m d (T) = 0.0134T- 3/11 m B H /n 



(34) 



We can then calculate the probability for a BH-WD sys- 
tem to contribute to the upper-end XLF at an ellipti- 
cal age T. We assume that (i) all accreting BH-WD 
systems were formed within a short interval of elliptical 
ages T s tart to Tjg n (in Gyrs) several Gyrs ago when star 
formation was still occurring in the elliptical; and (ii) 
T — Tfi n > i pcrs ,i.e., the binary is a transient at time T. 
The latter assumption is well justified given the short 
duration of the persistent phase (see §2.4). We further 
adopt a constant BH-WD formation rate between T sta rt 
and Tfi n , i.e., — const. The probability then is ex- 
pressed by the duty-cycle weighting factor at T for a 



W(T;m BH )= >..,, a* ^ 



f">d2 itt d 

. Q -^ JlTldl ™crit 



dt dm d u 



m^ m _Bt_ d 



Jmdi dt dm d 

CdT^d -1 * 7 ^ 



(35) 



Here m d i = m d (T - T start ; m BB ) = m d (ti; m BB ) and 
m d2 = m d (T - T Rn ;m BB ) = rn d (t 2 ,m BB ); m d2 > m dl . 
Then, using eq. l(T7|) . (|2T)|) and (|54*jl . we obtain: 



W(T; to B h) = 1-6 x lO^m^h 7/4 

1 - feAi)- 3/4 
l - (ta/ti) 



(36) 



Therefore for WD donors 7 — 0.5 and W depends on 
both the BH mass and the galaxy age. 

4. ACCRETING BLACK HOLE MASS SPECTRUM 

So far we have derived the dependence of the duty- 
cycle weighting factor W fea l2(j|) on the accreting BH 
mass and the age of the host elliptical galaxy, for the dif- 
ferent types of BH donors. In order to make progress and 
develop a method for deriving constraints on the slope (3 
of the underlying accreting BH mass spectrum we need to 
examine which of the possible donor populations domi- 
nate the observed upper-end XLF under what conditions. 
The answer to this question requires large-scale popula- 
tion synthesis models that are not part of the scope of 
this paper. However, as is shown below, we can use a 
number of different arguments and pieces of evidence to 
derive tentative constraints. The primary purpose of our 
analysis is not to derive an unambiguous constraint at 
present, but instead to develop a methodology for how 
to derive the most reliable constraints, given the current 
uncertainties associated with our current understanding 
of these X-ray binaries. 

For RG donors we find that in the case of a MT- 
dependent duty cycle 77 is about an order of magnitude 
smaller than for MS donors: by comparing eq. (|31|l and 
eq. (|33|l . using an elliptical-galaxy age of at least a few 
Gyr (more typical is 12 Gyr) and a BH mass of at least 
3M© the difference with the 77 value for MS donors is 
a factor of 15. We also note that among known BH X- 
ray transients in our Galaxy, the ratio of transient BH- 
RG systems to t ransient BH-MS systems is ab out 1:2 
(see Table 4.1 in IMcClintock fc Remillardl 12005]) . Con- 
sequently we conclude that BH-RG transients cannot be 
important contributors to the upper-end XLF of elliptical 
galaxies for the example case of the MT-dependent duty 
cycle as defined by 77. Only in the case of the constant, 
MT-independent (and therefore donor and BH-mass in- 
dependent) duty cycle we expect BH-RG transients to be 
a significant population of the observed upper-end XLF. 

Let us consider the case when the number of transient 
BH-RG systems exceeds the number of transient BH-MS 
system in a way that the contributions of the two popu- 
lations become comparable at some age of the elliptical. 
We also assume that /?o should be the same for both pop- 
ulations. In this case, the resultant combined XLF will 
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be flatter than the XLF provided by only BH-MS con- 
tributors. Secondly, as the contribution of BH-RG sys- 
tem decreases with elliptical age (see ea. l3*5|) . the XLF 
becomes steeper, evolving towards the slope character- 
istic for BH-MS binaries. It is possible that this is the 
kind of behavior that we observe in XLFs of ellipticals 
(see Fig. 3): we note that younger ellipticals appear to 
have flatter XLFs, although uncertainties are significant. 

For WD donors we find W{m) oc t^ 7 ^ 4 (see eq. OEJ), 
implying that the probability of each BH-WD transient 
contributing to the observed XLF decreases with the age 
of the elliptical galaxy (similar to the case of RG donors, 
but unlike the case of MS donors). Let us consider an 
elliptical where the formation of BH-WD systems has 
ended at least a few Gyr ago. We also consider that 
BH-MS systems are transient and have W(M) accord- 
ing eq. 1|31|) . In order for BH-WD binaries to contribute 
significantly to the observed XLF they must form at a 
rate such that more than ~ 8, 000 BH-WDs form for each 
BH-MS. This ratio is calculated adopting a value for the 
age of the elliptical of ~ 3.5 Gyr (among the lowest re- 
ported in the literature) and for a choice of BH masses 
in binaries with WD and MS donors, so that their Ed- 
dington X-ray luminosities are comparable, and there- 
fore they contribute to the same X-ray luminosity bin 
(3M Q for WD and 5M for MS donors). For more typ- 
ical, older ellipticals with ages closer to 10 Gyr the re- 
quired ratio becomes even higher than 8, 000. According 
to binary population synthesis models for the Milky Way 
published so far, the number of formed BH-WD LMXBs 
exceeds the number of BH-MS LMXBs by at most a fac- 
tor of 100 l|Hurlev et alJl2002D . Furthermore, the life- 
time of BH-MS binaries is of order 1 Gyr (or a few Gyr; 
see also Fig. 1), whereas the lifetime BH-WD binaries is 
longer, but cannot exceed the age of the elliptical galaxy 
(~ 10 Gyr). We conclude that the number of BH-WD 
LMXBs could be at most about a factor of 1000 higher 
than the number of BH-MS LMXBs, but this ratio is 
still below what is required for BH-WD to become an 
important contributor. So, if the ratio of BH-MS bina- 
ries to BH-RG binaries in ellipticals is similar to that in 
the Milky Way, BH-WD XRBs will not be a significant 
contributor to the XLF. Based again on the discrepancy 
between the duty cycles for WD and RG donors (smaller 
for WDs by a factor of ~ 800) , expect that BH-RG tran- 
sients dominate over BH-WD transients too. 

For the case of an example MT-dependent duty cycle 
(expressed by r\ in eq. |2J we conclude that: (i) if IT 
MB describes the angular momentum loss best, then only 
BH-MS transients significantly contribute to the XLFs 
of elliptical galaxies; consequently (3 = 2.5 ± 0.6; (ii) if 
instead RVJ MB is a better prescription, then the XLF 
is dominated by BH-RG binaries and (3 = 2.8 ± 0.6. 

For the case of a constant duty cycle independent of 
the donor type, it is clear that the XRB type with the 
highest formation rate should dominate t he XLF. Ac- 
cordin g to formation rates calculated by IHurlev et alJ 
(2002), BH-WD binaries form at a rate about 100 times 
higher than BH-MS and BH-RG binaries. Consequently, 
WDs would be expected to dominate the transient pop- 
ulation and this is certainly not true for the Milky Way. 
Therefore we conclude that the assumption of a constant, 
MT-independent duty cycle is most probably not realis- 



tic. 

Overall, we conclude that MS or RG donors dominate, 
depending on whether the IT or RVJ MB prescription is 
more realistic. Consequently, the slope of the accreting 
BH mass spectrum is ft = 2.5 ± 0.6 {f3 Q = 2.3 ± 0.6) or 
(3 = 2.8 ± 0.6, respectively. These quantitative results 
are of course dependent on the adopted example form 
of the MT-dependent duty cycle (eq. 0). As shown in 
§3.1.1 a possible different form (e.g., eq. PJ) could lead 
to somewhat flatter values for (3. 

Next we consider the fact that the upper-end XLF of 
ellipticals is not a perfect power-law up to arbitrarily 
high Lx values; instead there is a usually smooth cut- 
off behavior that limits the maximum Lx observed at 
~ 2 x 10 39 ergs -1 . In Fig.[3]we show the cumulative XLF 
associated with a model population of BH-MS binaries 
with a BH mass spectrum with a differential slope of 
(3 = 2.3 and with an imposed upper limit of 15 M Q on 
the maximum BH mass present in the XRB population. 
We obtain a model XLF that behaves very similarly to 
observed XLFs (dash-dotted line). Clearly this is just 
to show the importance of the qualitative effect of a BH 
mass cut-off on the cumulative XLF. 

5. DISCUSSION 

We consider the upper-end XLF of ellipticals (above 
the reported break at ~ 4 — 6 x 10 38 ergs -1 ) and suggest 
that it is populated by BH X-ray transients at outburst 
emitting approximately at the Eddington limit. We ar- 
gue that the upper-end XLF slope is a footprint of the 
underlying accreting BH mass spectrum modified by a 
weighting function related to the transient duty cycle. 
We show that this weighting factor is generally depen- 
dent on the BH mass and/or the age of the host galaxy 
and the derived power-law dependence is different for 
each of the possible BH donor types: MS, RG, and WD. 
Our predicted dominance of X-ray transients at outburst 
contributing to the upper-end XLF could possibly be 
tested by future high-resolution X-ray observations de- 
signed to achieve long-term monitoring probably at time 
scales of years or longer. Unfortunately, given the uncer- 
tainties in the theory of the thermal disk instability, it is 
not possible to make any predictions about the expected 
duration of these outbursts. 

Based on our analysis and prior population synthesis 
results we conclude that a constant transient duty cycle 
independent of the donor type can be excluded. Instead 
a duty cycle dependent on the binary and MT proper- 
ties seems to be required. Given the uncertainties asso- 
ciated with transient duty cycles at present, we adopt 
a couple of different formulations of such a dependence 
(eqs. |2I3| L as reasonable examples, which in no way ex- 
haust the possibilities. In the specific case of a duty cycle 
that depends on the ratio of the binary mass transfer rate 
to the critical rate for transient behavior (see eq. EJ we 
conclude find that the BH X-ray transients forming the 
upper-end XLF in ellipticals have a dominant donor type 
and an accreting BH mass spectrum slope [3 that depend 
on the strength of MB angular momentum loss: (i) for 
the IT MB prescription, only BH-MS transients signifi- 
cantly contribute to the upper-end XLF and (3 = 2.5±0.6 
(/?o = 2.3 ± 0.6 ); (ii) for the RVJ MB prescription, the 
XLF is dominated by BH-RG binaries and (3 = 2.8 ±0.6. 
We note that these quantitative results do depend on 
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our conclusions about which donor-type population dom- 
inates b ased on currently p ublished population synthesis 
models ( Hurley et al.l200"^) and on available observations 
of BH X-ray systems in our Galaxy. If, e.g., the relative 
fraction of BH-RG transients in ellipticals is larger than 
the observed relative fraction in our Galaxy, we expect 
that BH-RG binaries contribution will lead to a time- 
dependence of XLF slopes, where younger ellipticals will 
have a slope predicted for BH-RG binaries, and older 
ellipticals a steeper slop predicted for BH-MS binaries. 

The primary goals of this study are to present (i) the 
line of arguments that connects the upper-end XLF of el- 
lipticals to BH XRBs formed in the galactic field and (ii) 
the methodology for how to extract information about 
the accreting BH mass spectrum from the observed XLF 
slopes. We have further obtained quantitative results 
on the BH mass spectrum slope under certain reason- 
able assumptions, some of which (e.g., the functional 
form of the MT-dependent duty cycle) represent mere 
examples. A careful examination of the robustness of 
these quantitative results has been for the case of MS 
donors. It has been found that the derived slopes are 
robust against (i) random variations by factors of a few 
of the outburst peak luminosity of individual sources, 
and against (ii) variations of the possible duty-cycle de- 
pendence on the critical MT rate for outburst behavior. 
However, completely difference duty-cycle dependencies 
cannot be excluded. An improved understanding of this 



issue would be required to derive reliable quantitative 
conclusions about the value of the BH mass spectrum 
slope in transient XRBs in ellipticals. 

We expect that our analysis and methods can be used 
to reveal more information about the formation of BH 
XRBs in elliptical galaxies. More specifically they could 
eventually be used to constrain the physical connection 
between massive stars in XRB progenitors and the re- 
sultant BH masses. Current simulations assume either 
an artificially constant mass for BHs formed (usually at 
~ 10 Mq), or a constant mass fraction of the progeni- 
tor leading to the remnant objects, or a remnant mass 
relation consistent with core-collapse simulations. Con- 
straints on the accreting BH mass spectrum as those dis- 
cussed here could contribute to our understanding of core 
collapse, and the connection of BH masses to their pro- 
genitor masses. 
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